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Abstract: Two or more quantum systems are said to be in an entangled or non- 
factorisable state if their joint (supposedly pure) wave-function is not expressible as a 
product of individual wave functions but is instead a superposition of product states. 
It is only when the systems are in a factorisable state that they can be considered to 
be separated (in the sense of Bell). We show that whenever two quantum systems 
interact with each other, it is impossible that all factorisable states remain factoris- 
able during the interaction unless the full Hamiltonian does not couple these systems 
so to say unless they do not really interact. We also present certain conditions un- 
der which particular factorisable states remain factorisable although they represent a 
bipartite system whose components mutually interact and pay a particular attention 
to the case where the two particles interact mutually through an action at a distance 
in the three dimensional space. 



Introduction 



The term entanglement was first introduced by Schroedinger who described this 
as the characteristic trait of quantum mechanics, "the one that enforces its entire 
departure from classical lines of thought" [1]. Bell's inequalities [2] show that when 
two systems are prepared in an entangled state, the knowledge of the whole cannot 
be reduced to the knowledge of the parts, and that to some extent the systems lose 
their individuality. It is only when their joint wave- function is factorisable that they 
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are separable^. It is therefore interesting to investigate which are the situations such 
that two systems, initially prepared in a (pure) product state remain in such a state 
although they mutually interact. 

We shall show (sections 1 and 2) that when the Hilbert spaces associated to 
the interacting systems A and B are finite dimensional, if we impose that all the 
product states remain product states during the interaction, the full Hamiltonian 
can be factorised as follows: HAB{t) = HA{t) ®Ib + Ia® HB{t), where Hi{t) acts on 
the "-i" system only while Ij is the identity operator on the "j" system = A, B). 
In other words, in quantum mechanics there is no interaction without entanglement. 
We shall also present a sufficient condition under which particular factorisable (non- 
necessary pure) states remain factorisable during the interaction. 

We shall discuss, in the section 3, the situation where the interacting systems are 
two three-dimensional material points that interact through an action at a distance. 
We shall show that the factorisability of the full wave-function is preserved (A) in 
the test-particle limit (when one of the particles is quite more massive than the 
other one and is localised in a small region of space during the interaction), (B) 
in the classical limit which is considered here to play relatively to quantum wave 
mechanics a role comparable to the one played by geometrical optics relatively to 
classical wave optics and (C) when the Hartree approximation is valid. 



1 Two interacting spin one-halve particles 

Firstly, let us consider the most simple situation: the systems A and B are spin 
one-halve particles. We shall now show the following theorem: 

Theorem 0: 

Let us consider a system that consists of two spin one-halve particles A and B. 
Let us assume that the wave-function of the full system is a pure state of ® 
which evolves according to Schroedinger's equation: 

lhdt'^AB{t)=HAB{t)^AB{t) (1) 

^It can be shown that whenever two distant systems are in an entangled (pure) state, there 
exist well-chosen observables such that the associated correlations do not admit a local realist 
explanation, which is revealed by the violation of well-chosen Bell's inequalities [3,4]. 



2 



Let us assume that any arbitrary initially factorisable state ^j^si^ = 0) = ipA{t = 
0)^ipB{t = 0) remains factorisable during its temporal evolution: ^'as(^) = 
^ij(t) Wt > 0. 

Then, for each time t > 0, there exists a "factorisable" Hamiltonian H-^^^ = 
Ha ® /b + /a ® Hb (where Ha(b) is a self-adjoint operator that acts on C^) which 
brings the same change at time t as the change due to HABit)- 

Proof of the Theorem 0: 

Let us firstly consider that initially the state of the system belongs to a fac- 
torisable basis of ® that consists of the four following states: = 0) = 

\+a)®\+'b ), ^Asit = 0) = \+a) ®\-'b), ^Asii = 0) = \-a) ®\+'b), 
= 0) = I -a ) <S) I -'u ), where | + (')a(b)) and | - (')a(b)) represent up 
and down spin states along conventional axes of reference Za{Z'q) assigned to the 
systems A{B). 

Necessarily, ^\^{t) = \+a) ® I+b) for some couple of directions Za, Z'g. The 
unitarity of the evolution law imposes that either ^\^{t) = IV'K'^)) ® \~'b) 

= \~a) ® li^Bi't)) where \il'A{B)i^)) undetermined at this level of the proof. 
Let us consider now the first alternative and assume that the system is initially 
prepared in the product state ^ • + = I +a ) ® ^75 • (| +b ) + | )) . 

In virtue of the linearity of the evolution law, this state becomes at time t the state 
:^-(|+A)®|+'s) + |'i/'i(t))cg)|— '5)) which is not a product state unless = \+a) 

up to a global phase-factor that we can consistently take to be equal to unity (up 
to a redefinition of the phase of \—'b))- Then, ^ab{^) = 1+^) ® \~'b)- ^ similar 
reasoning, the second alternative leads to the conclusion that ^^^(i) = |— a)® 
By repeating this proof with "^ab instead of ^'^s, we get that either ^'^^(f) = 
|— a) ® I+a) or ^yiB(^) — I+a) ® \—'a)- virtue of the unitarity of the evolution 
law, must be orthogonal to ^'^^(t) so that, in conclusion, two alternatives 

remain possible: either (i) = I+a) ® p b) ^ab{^) = \~a) ® \+'b) or (ii) 

*ab(*) — \~a) ® \+'b) ^ab{^) — I+a) ® \—'b)- In ^-ny case, unitarity imposes 
that — \~a) (H) \ —'b) np to a global phase. If moreover we require that when 

the system is initially prepared in the product state |-(^'\^(0) + ^'^^(0) + ^'^g(0) + 
*ab(0)) it remains in a product state at time t, this global phase factor must be 
equal to unity. It is easy to check that in both cases all states that are initially 
product states are still product states at time t. 

Let us consider firstly that the first alternative is valid. It is easy to find a 
"factorisable" Hamiltonian H^/^\t) = HAit) ®Ib + Ia® HBit) that sends *^b(0) 
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on ^'^^(t) (i = 1,2,3,4) in a time t. Obviously, it is sufficient to choose Hji^(^B)if) 
in such a way that | +a{b) ) is sent onto | + in a time t. For instance we 

could take HA(B){t) to be a multiple of a time independent linear combination of 
the Pauli matrices that generates on the Bloch sphere a rotation that brings | +a{b) ) 
onto \+{')a{b)), and modulate the intensity of this Hamiltonian in order to perform 
the rotation in a time t. It is worth noting that such an Hamiltonian sends states 
that were initially product states on product states for all intermediate times t' 
{0 < t' < t). Note that we could even let depend Ha and Hb on time and tailor 
them in an ad-hoc way in order to generate arbitrary continuous state evolutions for 
all intermediate times t' (0 < t' < t) in so far the projections of \+Ait')) and 
are sufficiently regular curves that can be arbitrarily well approximated by a series 
of arcs of circles on the Bloch sphere. However, in the present approach, time was 
discretised, and it is not so simple to quantify properly what we mean by "regular" , 
a limitation that we shall overcome in the next section, in the proof of our main 
theorem, where the hypothesis of regularity in time is expressed quite naturally by 
a requirement on the Taylor development of the temporal state evolution. 

Let us now consider the second alternative and assume that for all intermediate 
times f (0 < t' < t) the evolution sends states that were initially product states 
on product states. Then, at time t' — ^ either the first alternative is valid or the 
second alternative is valid. 

If the first alternative is valid at time t' let us consider the time t" — ^. At time 
t" either the first alternative is valid or the second alternative is valid. If the first 
alternative is valid let us consider the time t'" = ^. Otherwise let us consider the 
time t'" = f . 

If the second alternative is valid at time t' let us consider the time t" = |. At 
time t" either the ffist alternative is valid or the second alternative is valid. If the 
ffist alternative is valid let us consider the time t'" — ^. Otherwise let us consider 
the time t'" = g and so on. 

By doing so it is easy to show that for any positive integer there must exist 
an intermediate time to such that the states ^^si^ = ^o) = I +a )o ® I +b )o' 

^Asit = to) = \+a)o^\ -b )o' ^ABit = ^o) = I -A )o ® I +'b )o' ^Asi^ = 

to) = \-a)o<S)\ -'b )o> (where | + {^)a{b))q and | - C)a(s))o represent up and down 
spin states along conventional axes of reference Z\[Zb) assigned to the systems 
A{B)) are sent at time to + e onto the states ^^^^(t = to + e) = \ +a ) ^ ® \ +'b )o 

"^Isit = to + e) = \-A)M+'B)e^ ^ABit = to + e) = \+A)M-'B)e^ = ^0 + ^) = 

I -A )e ® I -'b )e' with e = ^ (where | + (')A(B))g and | - (')A(s))g represent up and 
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down spin states along conventional axes of reference Z^(Z^). For sufficiently 
large, the continuity of the evolution law (eqn.|I]) imposes that | +A)e = I ~^a)o + t{^) 
and I = I +b)o + t{^) where by definition r(e™') decreases at least as fast as the 
mth power of e when e goes to zero. But then the in- product between ^^^(t = to) 
and ^^^(t = te) is equal to zero (up to T(e)). In virtue of Pythagoras's theorem 
their distance must be equal to \/2 (up to r(e)) which contradicts the continuity of 
the evolution law. Therefore, the second alternative is not valid and for any time t, 
we can find a "factorisable" Hamiltonian HABif) = HAif) ® Ib + Ia® Hsif) that 
sends ^^^(0) on ^^^{t) {i = 1,2,3,4) and sends product states on product states 
for all intermediate times. 

Note that continuity in time plays a crucial role in our proof which is similar to 
the proof of the impossibility of sending a righthand glove onto a lefthand glove by 
a continuous isometry of the three dimensional physical (Euclidean) space, an intu- 
itively obvious fact that presents deep analogies with the situation encountered here. 
The transformation described in the second alternative can effectively be obtained 
by composing the transformation described in the first alternative with a discrete 
transformation during which the states of the systems A and B are interchanged. 
Similarly, the composition of a continuous isometry (composition of Euclidean ro- 
tations and translations) and of a refiection in a mirror sends a righthand glove on 
a "virtual", non-physical lefthand glove. It is highly probable that we could refine 
the theorem and generalise it to arbitrary finite dimensional systems (where the 
full system is represented in the Hilbert space C^-* C^, with dA and ds standing 
for the dimensions of the systems A and B) by following the same way of reasoning. 
Nevertheless, we shall now give a general proof based on a different approach that 
is less abstract but appeals more to physical intuition. 



2 Two interacting finite-dimensional systems 



We shall now present our main result. Let us consider two interacting quantum 
systems A and B. We assume that the Hilbert spaces associated to these systems 
are finite dimensional (of dimensions rf^ and ds respectively), that the wave-function 
of the full system is a pure state of C^^ ® C'^^ and obeys the Schroedinger equation: 
ihdt^ABii) = HAB{t)'^ABit) where HABif) is a self-adjoint operator that acts on 
C'^-'^ ® C^^, that we assume to be sufficiently regular in time in order to ensure 
that the temporal Taylor development of the wave-function is valid up to the second 
order in time. 
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Main Theorem: 

All the product states remain product states during the interaction if and only 
if the full Hamiltonian can be factorised as follows: 

HAB{t) = HA{t) ®Ib + Ia® Hsit) (2) 

where acts on the zth system only while Ij is the identity operator on the jih. 
system (i, j — A^B). 

In order to prove this theorem, we shall firstly prove the following lemma: 

Lemma: 

A pure product state remains product state during the interaction if and only if, 
during its evolution, the Hamiltonian never couples this product state to a product 
state that is bi-orthogonal to it. 

Proof of the Lemma: 

A) Proof of the necessary condition. Let us consider that at time t the system 
is prepared in a product state ^'ab(^) = i^Aif) ®i/^B{t), and let us choose a basis of 
product states ® {i : l...dA;j ■ and {i/^aI'^a) = % = (^sI^b)) such 

that ipA{t) — iV'i) ^-nd (t)B{t) — Then, after a short time 5t, 

iSf 

+ 5t)^{I+—- HAB{t)) ■ ^AB{t) + T{8t^) 

where r(e™) was defined in the previous section. In a matricial form, the previous 
equation becomes: 

i5t 

where 

H^kJl = m ® {<P^\HAB{tM'A) ® M + r{Sf) (4) 

Equivalently, 

iSt 

^AB{t + 6t) = ® 10^) + —{^rA...d^Hnil\i^\) ® 10^) 
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All the components of ^ AB{t + St) that are bi-orthogonal to ^Asit) are contained in 
the last term of the previous equation: Sj;2...dA;i:2...rfs-^jiii|'0A) ® l^s))' ^^P to r(5f^). 
We can rewrite this equation as follows: 

iSt i6t 

i6t 

+ -J-^^■.2...dA■,r■2...dnHil,l\^p\) ® + r(5t2) (5) 

Let us assume that the Hamiltonian couples ^AB{t = 0) to states that are bi- 

2 

orthogonal to it, which means that ^i:2...dA;j:2...dB 

\Haj^\^ ^ 0. We shall now show 
that then the development of the first order in 6t of the bi-orthogonal or Schmidt 
decomposition [5] of '^Asi't + St) contains more than one product state, which means 
that ^ AB{t+St) is entangled for 5t small enough. In order to do so, let us consider the 
components of ^Asit + St) that are bi-orthogonal to ^Asit)- In virtue of Schmidt's 
theorem of the bi-orthogonal decomposition [5] , we can find dA — l normalized states 
\iPa) of C*^^ mutually orthogonal and orthogonal to \iPa) and ds — ^ normalized states 
\(f)^'^) of C"'^ mutually orthogonal and orthogonal to \(f)^^) such that 

i6t ■ ■ 

-^^i:2...dA;j:2...dBHiljl\lp\) ® = Si:2...m*n(d^,dB)"#A) ® 

Let US now define the state of C^-^ as follows: \ip]l) = ^■(|V'A) + x^»:i-'^A-^iiii|^A)), 
where A'"! is a normalisation factor, and let us replace the orthonormal basis iV^i), 1"^!), 

of C^^ by the orthonormal basis {\'i/ja) ,\'4'a) A'^a) 1 ■■■^ I'^l*)} of ^'^^ that we obtain 
by the Gram-Schmidt orthonormalisation procedure: 

lv^i') = ^-(l^i)-(^Al^i)-IV'i')) 

where A'2 is a normalisation factor. 

I^i') = ^ ■ (l^i) - (^i'l^i) ■ \^a) - (V^l'l^i) ■ \^a)) 

where is a normalisation factor, and so on. It is easy to check that \iPa) = 
+ T{5t). Note that this is no longer true when the dimension dA is not fi- 
nite. We can repeat the same operation in order to replace the orthonormal basis 
Wh), 1^^^)} of C^- by the orthonormal basis |^|), |^,'^)} 

of C^^. Then, after substitition in the eqn.^, we obtain that: 

*AB(t + 5t) = ^i:l...min{dA,dB)C^i\i^A) ® \<Pb) + '^i^t^) 
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where 



The previous equation expresses that the development up to the first order in 6t 
of the bi-orthogonal Schmidt decomposition of ^AB{t + St) contains more than 
one product state. It is well known that then ^Asit + St) is an entangled state. 
Nevertheless, for those who are not familiar with this property, we shall prove 
directly the result by making use of the reduced density matrix. By definition, 
the reduced density matrix of the system A is equal to Tr^p where p is the 
projector on ^ab- Obviously, when the state of the system is a product state 
{^AB = i^A ^ i^B), Pa is the projector on i/ja, and we have that pA = Pa, and 
TrpA = Trp\ = 1. Actually, TrpA - Trp\ provides a good measure of the degree of 
the entanglement of the full system. If the Schmidt bi-orthogonal decomposition of 
the state ^ ab is equal to '^i:i...min{dA,dB)^'i\''P^A) ® Wb)i then it is easy to check that 
Pa = T,iA...rninidA,dB)\'^if\^A){'^A\i TvpA = ^i:i...min{dA,dB)\'^if = 1 by normalisation 
and Trp\ = ^i:i...min{dA,dB)\'^i\'^ — (TrpA)"^ = 1^ = 1. The last inequality is saturated 
for product states only. Note that Trp\ = Trp\ which shows that this parameter 
expresses properties of the system considered as a whole, as it must be when we are 
dealing with entanglement. Obviously Trp\{t + 5t) < \ai\^ + {X'i:2...min{dA,dB)\'^i\'Y ■ 

But \ai\^ = (1 - ^i:2...'min{dA,dB)Wi?f = (1 - ^Si:2...d^;j :2...ds | i^iljl | ^ + T^dt^))"^ and 

{'>^i:2...min{dA4B)Wi\^f = i^^i:2...dA;j:2...dB\Hiijif + r{5t^))'^ in virtuc of the eqn.| so 
that Trp\{t + St) < 1 - 2 ■ ^J:i.,2...dA;j:2...dB\Hiijif + T{6t^) < 1 for St small enough^, 
which proves the necessary condition of the lemma. 

B) Proof of the sufficient condition. Let us consider that at time t the system 
is prepared in a product state ^'ab(^) = V'a(^) ® i^Bit), and let us choose a basis 
of product states I^Pa) ® I^b) similar to the basis introduced in the proof of the 
necessary condition. When the Hamiltonian does not couple ^AB{t) to states that 
are bi-orthogonal to it, ^i:2...dA;j-2...dB\Hiijif = (where Hikji is defined in the eqn.|) 
and, in virtue of the eqn.|^: 

ihdt^ffABit) = HAB{t)^AB{t) = (S,:i...,,i/,ni|V'A))®l'^B) + l^i)®(Si:2...dsi^ll,l|<^i)) 

(7) 

■^It can be shown by direct computation that when the state of the system is a product state 

{^Asit) = ^Ait) (8) ipB{t)), then the foUowing identity ^^^^(0 = is necessarily satisfied, inde- 
pendently of the form of the Hamiltonian Hab- This explains why no term of the first order in 
time appears in the previous development. 
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We can rewrite this equation as follows: 

ihdt ^ABit) = m"-it) ■ Mt)) ® Mt) + Mt) ® (HB^it) ■ Mt)) (8) 

where the effective Hamiltonians H'^-^-l^' are defined as follows: 

H'J^it) ■ pA{t) = TTB{HAB{t)pAB{t)) (9) 

and 

H'li^-it) ■ pB{t) = TTA{HAB{t)pAB{t)) - {TTAB{HAB{t)pAB{t))) ■ PB{t) (10) 

In these expressions Tvi represents the partial trace over the degrees of freedom as- 
signed to the system i while pAB{t) is the projector onto Pa(^) = TrBPAB{t), 
and pB{t) — TrAPAB^t)- For instance, we have that 

TrB{HABit)pAB{t)) = ^,.i...dMHAB\^\) ® ® 

= Si:i...<iB(0B|Si:l...<iAj:l...<iB^^'ajl|'0A) ® M^nM 
— ^l:l..MB'^i:l---dA,j:i---dBHiljl\'lpA)^lj^ll{'^A\ 

= (S,:i...rf,i7aii|V'i>(V'\l) 
so that H'^"-{t) ■ iljA{t) = ^i:i...dAHiiiM- 

Let us consider the product state ipA^^it') (8) ip^'^{t'), where ipA{B){t') ^ ^O" 
lution of the reduced Schroedinger equation ih dt' i/j^^^^it') — H^/f^g^{t') • V'a(b)(^') 
for the initial condition = i-' a(B){^) ■ Obviously, ih dt' 'ip'^^\t') ® tp'^'^^{t') = 

HABVf{t') ® VB^t') and ^ AB{t) = ij^'l'^it) (g) ijj's^^{t) so that *ab(^') = i^T'^it') ® 
ip^B'^{t'), Vt' > t which ends the proof of the lemma. 

We shall now prove the main theorem. 

Proof of the Main Theorem: 

A) Proof of the necessary condition. Let us choose a basis of product states 

li^A.) ® {(Pb) • ^■■■dA',j '■ l---dB and (V'aIV'a) = ^ij = {'^bI'^b))- If impose that 
all the product states remain product states during the interaction, then, in virtue 
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of tlie lemma, the full Hamiltonian never couples a product state to a product state 
that is bi-orthogonal to it. Then, at any time t, ^i:2...dA;j-2...dB\Hiijif = (where 
Hikji is defined in the eqn.^ so that we have that: 

where 

|A^M) = S,.i,..,,i7,,,,#l) (11) 

and 

I^B'Pb) = ^k:l...dB,k=^jHiikj\(j)B) (12) 

Let us consider that at time t the system is prepared along one of the first four 

states : 1, -4) of this basis: ^^{t) = I^a) ® \<Pb), = I^A) ® \<f^l), 

^\Bit) = m ® 101,), ^ut) = m ® m- Then, 

HMt) ■ = ® \^b) + \^\) ® I^b^b) 

//^^(t) ■ *i^(t) = lA^Vi) ® M + l^i) ® |A|V1,) 

HABit) ■ *%(t) = lAiVi) ® 102b) + l^i) ® |A|20|) 

By linearity, 

4((|A^Vi) + |A^V^)) ® 10^) + ® |A1,V^) + m ® |A|V^)) 



1 

V2 



4((|A^Vi) + lA^Vl)) ® 10^) + (l^i) + 1^1)) ® (|A^V^) + lAlVl,)) 



+ (lV^i)-|^^))®(|A^V^)-|A|V^)) 

TldV^i) - iV^i)) is orthogonal to ^(|^i) + |^i)), so that HABit) ■ j^.i^Uit) + 

*ab(0) couples -^i^ABit) + ^ABit)) to a bi-orthogonal state unless (|A^Vb) — 

I A^^0b)) is parallel to \4>b)- Now, :^(^ab(0 + ^ABit)) is a product state so that, 
in virtue of the lemma, the following constraint must be satisfied: 

H^'bVb) - I^IVb)) = A|0^) 
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The same reasoning is valid with the states ■:^{^ABi^) + ^abI^))? ':^i^ABi'^) + 
*a_b(^)) '^i^ABi^) + ^AB^'^)) leads to the following constraints: 

(|A]|0|) - \A^i<Pl)) = AVI) 



(lAiVi) - |A-^1)) = A-IV'^) 

By definition (eqn.|l2]), | A^0;^) is orthogonal to |(;/)'^) so that necessarily A = A' = 
0. Let us now consider the product state (*ab(0 + ^ABi^) + ^ABi^) + 
By linearity: 

HAB{tyl{^Ut)+^Ut)+^Ut)+^Ut)) = 

= -^((|A\Vi) + lA^VD) ® 101) + (|AM) + lA^^^i)) ® 101) 

+ ® (I A|VI) + |A|20|)) + 1^^) ® (|A|VI) + |A|VI))) 
In virtue of the constraints, we get that: 

= + ^"1V4)) ® \<Pl) + i\^Ai^\) + I AM)) ® (I'/'b) + 14)) 

+(l^i) + l^i))®(|A|Vl.) + |A]|0|))) 

Such a state does not contain any state bi-orthogonal to + ^ab(0 + 

^^^(t) + ^Ui^)) only if A"|7/;i) + = X""{\^\) + |^i)), which imposes that 

A" = A'" = A"". We can repeat this proof with the indices ii' for the system A and 
Ij for the system B instead of 12 as it was the case in the previous proof, and we 
obtain that |A^i4) = |A^;^0i) = jA^^i), and jA^^'^^) = |A^Va) - A^OI^ji) = 
|A^'0^) — A(j)|'0^) (where IA^t/^^) does not depend on j while A(j) and |A-^(/)-^) 
do not depend on i). Therefore: 

HABit) ■ IV'a) ® \<l>i) = l^Ai^A) ® I0'b) + I^a) ® |Ai,0i,) - A(j)IV'a) ® I0'b) 
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which fulfills the eqn.^ provided we proceed to the following identifications: H^lt) ■ 
= \Ai^^/j\) and //^(t) ■ {(j)^) = \A{(p{) - X{j)\(j)i). This ends the proof of the 
necessary condition of the main theorem. 



B) Proof of the sufficient condition. Let us assume that the full Hamilto- 
nian can be factorised according to the eqn.||. Let us consider the product state 
ipA'^it') where ipA{B)i'^') a solution of the reduced Schroedinger equation 

ihdt> ^7(13) (f) = i/A(B)(t')-V^Afi?)(^') for the initial condition VA{B)ii) = '^A(B)it). Ob- 
viously, ihdt'ip''f{t')®^]^'^{t') = HAsit'WA'^it') ^ilj's^'^{t') and = Vl'^it) ® 
V'B'^(t) so that ^Asit') = i^A'^it') ® Vs^it'), \/t' > t which ends the proof of the 
main theorem. 



Note that the condition ^ encountered in the proof of the sufficient condition 
of the lemma can be generalised to factorisable non-necessarily pure states. This is 
the essence of the following theorem. 

Theorem 2: 

If initially, a bipartite system is prepared in a factorisable (non-necessarily pure) 
state: pAsit = 0) = pAit = 0) ® psii = 0), and that Vt > 

HAB{t)pAB{t) = {H'J^it) ■ pA{t)) ® pB{t) + PA{t) ® {H's^^-{t) ■ pB{t)) (13) 

where 

H'J^it) ■ PAit) = TrB{HAB{t)pAB{t)) 

and 

H'J^-{t) ■ pBit) = TrA{HAB{t)pAB{t)) - {TrAB{HAB{t)pAB{t))) ■ pB{t), 

then, necessarily, the state remains factorisable during the interaction: pABit) = 

PAit)® PBit) Vt> 0. 

Proof of the Theorem 2: 

When we describe the state of the system by a density matrix, its evolution 
obeys the von Neumann equation: 



in dt PABit) = [HABit), PABit)] (14) 

where [^,5^] represents the commutator of two operators X and Y. HABit) and 
PABit) are self-adjoint operators so that [HABit), PABit)] = HABit)- PABit) — iHABit)- 
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PAB{t))~^ where O"*" represents the self-adjoint operator of O. Moreover, in virtue of 
the eqn.|T3|, we have that: 

ihdt pABit) = {H'^J^-it) ■ paH)) ® psit) + PaH) ® (H'J^-it) ■ psit)) 

- HHff-it) ■ pAit)) ® psit) + pAit) ® {H'^JHt) ■ PB{t))y (15) 

Let us consider the product state p^X'^if) ® p^'^(t), where p^'^(B)(t) is a solution 

of the reduced von Neumann equation ih dt p'^X(B)i^) ~ ['^A(B)i^)^PA{B)i^)] 
initial condition p^fB)(t = 0) = pA(B){t = 0). Obviously, ihdt p''X'^{t) ® Pe'^it) = 
HAB{t)pTit) ® Ps'^it) and pAsit = 0) = pj''{t = 0) ® = 0) so that pAsit) = 

PA'^i't) p^s'^it), Vt > which proves the theorem^. 



Note that the eqn.|T3| is linear in the coupling Hamiltonian Hab and is auto- 
matically satisfied when the eqn.^ is satisfied. Nevertheless it is non- linear in pab- 
Moreover, the effective potential that acts onto say the A particle is likely to depend 
on the state of the B particle, a situation that does not occur if we impose that all 
product states remain product states. 



Beside, it is worth noting that the sufficient condition expressed by the eqn.|T3| is 
also necessary in the case of pure states. Effectively, if, initially, the bipartite system 
is prepared in a factorisable pure state: ^AB^t = 0) = = 0) ® il^B{t = 0), 

and that the state remains factorisable during the interaction, then, in virtue of 
the necessary condition of the lemma, the Hamiltonian may not couple the state 
^AB{t) = i^Aif) at any time t > to a bi-orthogonal state so that, following 

the proof of the lemma, the eqn.^ must be valid at any time. Therefore, in virtue 
of the eqn.[^, the eqn.|l3| must be valid too. 



Now, the sufficient condition expressed by the eqn.[T^ is in general not necessary 
in the case of non-pure states as shows the following counterexample. If initially, the 
bipartite system is prepared in a factorisable state: pAB(t = 0) = pA^t = 0) = 
0), and that Wt > 0, HABif) = PAB^t = 0), then it is easy to check that pAB^t = 0) = 
PABim > 0, ■ PA{t) = TrB{HAB{t)pAB{t)) = p\{t = 0) ■ Tr^pK^ = 0), 

W • PB{t) = TrA^HABit) PABit)) - {TrAB{HAB{t) PABit))) ■ pnit) = TrAp\{t = 

0) . = 0) - Tr-ApAit = 0) ■ TrBp\{t = Q) ■ p^it = 0) and it is easy to check that 
in general the eqns. || or |l^ are not valid when the initial state is not pure so to say 
when it is not a product of pure states. 

''Note that this proof as well as the proof of the sufficient condition of the main theorem are 
also valid when the systems A and B are infinite dimensional, for instance when they are localised 
particles that interact through a central potential. 
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3 The infinite dimensional case 



The proofs of the necessary conditions of the lemma (and thus of the necessary 
condition of the main theorem) are not valid when the systems A and B are infinite 
dimensional. Nevertheless, we conjecture that these conditions are still true in that 
case, so to say that there is no interaction without entanglement. Let us for instance 
consider that A and B are two distinguishable particles, and that their interaction 
potential is an action a distance that is time-independent and invariant under spatial 
translations (a Coulombian interaction for instance). They fulfill thus (in the non- 
relativistic regime) the following Schroedinger equation: 

ihdt'^irA, vb, t) = -(- Aa + Ai3)^(r^, r^, t) 

Ztua ZniB 

+ VAB{rA - rB)*(rA, r^, t) (16) 

where A.a{b) is the Laplacian operator in the A{B) coordinates. As the potential 
does depend on the relative position r^ej = — only, it is convenient to pass to 
the center of mass coordinates: 

iUdt'^iYcM, Yreh t) = "(tTT \ tAcM + — A^e/)*(rcM, r^ei, t) 

2{mA + mB) 2fj, 

+ V^B(l"re/)^(l"CM5 T^rel, t) (17) 

where tcm = "'"'''"'T'"-^"'' and u = iiM^. As it is well-known, the previous 
equation is separable which means that if, initially, the wave-function is factoris- 
able in these coordinates, it will remain so during the evolution. Now, we are 
interested in sitTiations for which the wave-function is initially factorisable accord- 
ing to the partition of the Hilbert space that is induced by the systems A and 
B. In general, such a wave-function is not factorisable in the coordinates of the 
center of mass. Formally, if ^(r^, r^, t = 0) = •0a(i'a, t = Q) ■ 4>b{yb, t = 0), 
*(rcM, Yrei, t = 0) = JdLuA{Lu)'ilj'^j^^{rcM, t = 0) • il^reii^reh t ^ 0) where A{uj) is 
a generally non-peaked amplitude distribution. Then, at time t, "^{tcm, T^rei, t) — 
J dujA{u)ij'^Mi^CM, t) ■ ijj[;^^i{rrei, t) , whcrc ip'^j^,j{rcM, t) obeyed during the time in- 
terval [0,t] a free Schroedinger evolution for the initial condition 'il)'^j^{rcM, t = 0) 
while '0rei('^reZ) i) was submitted to the interaction potential VABi^rei)- In general, 
^(i"a, i"b, t) is no longer factorisable into a product of the form V^A(i"yi, t)-'4>B{^B, t). 
Actually, this is not astonishing because, in virtue of Noether's theorem the full mo- 
mentum is conserved during the evolution. Therefore the recoil of one of the two 
particles could be used in order to determine (up to the initial undeterminacy of 
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the centre of mass) what is the recoil of the second particle. The existence of such 
correlations is expressed by the entanglement of the full wave- function. On the basis 
of such general considerations we expect that entanglement is very likely to occur 
due to the interaction between the two particles. 

Nevertheless, it is interesting to investigate in which situations it is a good 
approximation to consider that the systems A and B remain in a factorisable state 
during time. We shall distinguish three typical situations. 

3.1 Scattering of a light particle by a heavy and well local- 
ized target (the test-particle limit) 

Let us assume that << m^, and that the B particle is initially at rest and well 
localized. The particle A is assumed to be initially prepared in such a way that it 
will pass in the vicinity of the heavy particle B, that its trajectory will undergo a 
deviation due to the influence of the interaction Vab, and that it will finally escape 
to infinity without exerting any significant back action onto the particle B. This 
situation is often encountered during scattering experiment. If we let coincide the 
origin of the system of coordinates associated to the particle B with its location, and 
that we neglect its recoil as well as its dispersion (this approximation is only valid 
during a limited period of time), the following approximations are valid: tcm ~ 

Tb ~ 0, Trel - = T^, '0A(rA, t) ^ Ipreli'^reh t) and V's(rB, t) ^ ^CAf(rCM, *)• 

Moreover, ^(r^, r^, t = 0) = ipA{rA, t = 0) ■ V'i?(rB, t = 0) ^ Aei{rrei, t = 0) ■ 

'4'CM{rcM, t = 0) ^ ^(rcM, r^ei, t = 0). At time t, ^'(rcM, ^rd, t) ~ Aeli^rel, t = 

0) • i^cMiycM, t) pa ilJAiyA, t) ■ ^As(rs, t) ^ "^{va, re, t). The separability of the full 
system into its components A and B is thus ensured, in good approximation, during 
the scattering process. 

3.2 Mutual scattering of two well localized wave packets 
(the classical limit-interacting material points) 

Another interesting limiting case is the situation during which we can neglect the 
quantum extension of the interacting particles. This will occur when the interaction 
potential Vab is smooth enough and that the particles A and B are described by 
wave packets the extension of which is small in comparison to the typical lenght of 
variation of the potential. It is well known that in this regime, when the de Broglie 



15 



wave lenghts of the wave packets are large enough, it is consistent to approximate 
quantum wave mechanics by its geometrical limit, which is classical mechanics. 
For instance the quantum differential cross sections converge in the limit of short 
wave-lenghts to the corresponding classical cross section. Ehrenfest's theorem also 
predicts that when we can neglect the quantum fluctuations, which is the case here, 
the average motions are nearly classical and provide a good approximation to the 
behaviour of the full wave-packet in so far we consider it as a material point. In 
this regime, we can in good approximation replace the interaction potential by the 
first order term of its Taylor development around the centers of the wave-packets 
associated to the particles A and B: 

VAB{rA-rB) ~ Vab{< >t - < >t)+VAVAB{< >t - < rs >t)-{rA- < ta >t) 

+^bVab{< ta > - < rs >t) • {tb- < tb >t)- 

Then the evolution equation is in good approximation separable into the coordinates 

r_4,rB and we have that, when \f(rA. r^j, t = 0) = ipAi^A, t — 0) ■ '4>b{'<^b, t — 0), 
then, at time t, ^(r^, r^, t) ~ t) ■ iPb{^b, t) where 

ihdtipA{rA, t) ^ AAipA{rA, t) 

Zuia 

+ {Vab{< rA >t> - <rB >t)+'^AVAB{< rA>t- < rs >t)irA- < ^a >t))'ipA{rA, t) 

(18) 

^2 

indti^BirB, t) fa ABipBirB, t) 

2mB 

+ {"^bVab{< VA>t- < Vb >t) ■ (rs- < Tb >t))i^B{rB, t) (19) 

Note that the Bohmian velocities associated to the particles A and B are factorisable 
only when the full state is factorisable. Otherwise, the velocity of a particle depends 
non-locally on the location of both particles. 



3.3 Bound states: the Hartree approximation 

When the energy of the full system is negative, we expect that it will remain in a 
well localised bound state. When one particle is quite more massive than the rest of 
the system as is the case with the sun in the solar system or with the nucleus inside 
the atom, it is a very good approximation to neglect its recoil and its extension 
for all times. Indeed, if we think to the nucleus for instance, its recoil is zero in 
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average over an orbit, and its Compton wave lenght is very small. Therefore it is 
consistent in a first approach to reduce the study of the energy levels of atoms to the 
study of the energy levels of the electrons that are assumed to undergo an external 
central Coulombian potential due to the presence of the nucleus and to factorize the 
full wave function into a product of an electronic wave function and of a nuclear 
one. For sure this approximation is valid to the extent that we can neglect other 
degrees of freedom as the nuclear spin and so on which is not always the case. If 
moreover we assume that the electronic wave function is itself factorisable, which 
is certainly a crude approximation because of the presence of exchange terms due 
the undistinguishability of the electrons and because the Coulombian interaction 
between the electrons is likely to generate entanglement, we are performing the so 
called Hartree| approximation [6]. Let us consider the Helium atom for instance, and 
let us neglect the fermionic exchange contributions, the spins of the electrons and 
of the nucleus and so on. The time independent (electronic) Schroedinger equation 
is then the following: 

EAB-'^irA, tb) = Aa + Va- Ab + l^B)^(rA, r^) + VAB{rA-rB)'^{rA, r^) 

ZTJIa ZTTLb 

(20) 

where Va and Vb represent the external fields (for instance the Coulombian nuclear 
field), while Vab represents the Coulombian repulsion between the electrons A and 
B. Let us assume that this equation admits a factorisable solution ^'(r^, r^) = 
V'A(rA) ■ i^B{rB)] then: 



Eab ■ ipA^rA) ■ ipB{rB) 

= {{-- Aa + VA)^A{rA)) ■ i^Bi^B) + ^A{rA) ■ {-- + VB)^{rB) 

ZrriA ZTfiB 

+ Vab{^a - rB)^A(rA) ■ V'B(ri?) (21) 

Let us now take the in-product of this equation with iPa{j^a) and multiply the 
resulting equation by ^1)a{ja) respectively. We obtain: 



Eab ■ ^A(rA) ■ i^Bij^B) 

^Note that when the Hartree approximation is vahd, particles behave as if they were discernable, 
and constituted of a dilute, continuous medium ditributed in space according to the quantum 
distribution in V'a(b)(''A(S)i which is close to the interpretation of the wave- function originally 
adopted by Schroedinger. 
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= ^^M- < {-- + Va) >a -i^BirB) + i^Ai^A) ■ (-7^ + VB)^|J{rB) 

ZUIa ZTTIb 

+ < VAsirA - rs) >A -ipAirA) ■ ipBirs) (22) 
Similarly, we get that: 

Eab ■ ^A(rA) ■ V'slri?) 
=< (-7^ Ab + Vb) >b -iJAirA) -i^BirB) + {{-- Aa + V:i)^(rA)) ■ V^slrs) 

ZTTIa ZUIb 

+ < VAB(rA - Vb) >b •^A(rA) • V^B(rB) (23) 



Let us now sum the two previous equations and substract the eqn.^. We obtain 
the following consistency condition: 

{Eab- < {-- Aa + Va) >a-< A^ + Vb) >b) ■ ^^(rA) ■ ^^(rij) 

ZniA ZniA 

= (< VABi^A-VB) >A + < VABi^A-VB) >B -V^Aij(rA-ri?)) ■^^(rA) ■'^^(rB) (24) 

Equivalently, when the wave-function does not vanish, the following condition must 

be satisfied: 

Vab{va - Vb) =< Vab{va - vb) >a + < Vab{va - vb) >b - < VAsi^A - vb) >ab 

(25) 



This is nothing else than the condition 13 in a static form. Note that here it appears 



to be a necessary condition , which does not infirm our conjecture that the necessary 
condition of the lemma is true in the infinite dimensional case. Reciprocally, it is easy 



to check that if the eqn.|T3| is satisfied, that the full state is factorisable and that the 
reduced states of the particles A and B are eigenstates of their respective effective 
Hamiltonians, the full state is eigenstate of the full Hamiltonian, in agreement with 
the sufficient condition of the lemma (which is a special case of the theorem 2). We 
see thus that the Hartree approximation is valid when the interaction factorises into 
the sum of two effective potentials that act separately on both particles, and express 
the average influence due to the presence of the other particle (which is not true 
in general and certainly not inside the atom). Similarly, in the test-particle limit, 
the effective potential undergone by the massive particle is close to zero, and when 
the heavy particle is well localised, its average, effective, potential is close to the 
real potential undergone by the light "test-particle" . In the classical limit (material 
points), the quantum internal structure of the interacting particles can be neglected 
and the potential is equivalent to the sum of the effective potentials. 

In the three cases, the necessary condition of the lemma is not infirmed, and its 
sufficient condition is confirmed. In all the cases, the systems are separated only in 
first approximation. 
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4 Conclusions and comments 



A conclusion of this work could be: in quantum mechanics to interact means nearly 
always to entangle. We showed that real interactions do necessarily generate entan- 
glement (the inverse result, that it is impossible to generate entanglement without 
turning on an interaction, is rather trivial). 

Considered so, the degree of entanglement of the universe ought to increase with 
time, which would indicate some analogy between entanglement and entropy. Note 
however that the temporal reversibility of the Schroedinger equation implies that the 
degree of entanglement could also decrease in time so that we face a paradox ana- 
log to the famous Loschmidt paradox which emphasises the apparent contradiction 
between the temporal asymmetry of the second principle of thermodynamics and 
the temporal symmetry of fundamental interactions. Obviously, such considerations 
are out of the scope of this paper and we invite the interested reader to consult the 
reference [7] and references therein. 

Beside, it would be worth investigating the generalisation of our results to infinite 
dimensional systems. We conjecture that the necessary conditions of our lemma and 
of our main theorem are still true when we deal with infinite dimensional systems, 
as is the case for the corresponding sufficient conditions and for the theorem 2. 

Let us briefly reconsider the three situations during which the interaction be- 
tween two mutually interacting particles (in the three dimensional, physical space) 
is entanglement-free, at least in first approximation (see section 3). These are the 
test-particle limit (no feedback), the geometrical limit of quantum wave mechanics 
(narrow wave-packets) and the Hartree approximation (particles seen as a dilute 
gas). Each of these situations has a counterpart in "classical" physics: idealised 
test-particles play an important role in classical mechanics and in general relativity, 
the geometrical limit of quantum wave mechanics is Hamiltonian mechanics, while 
the image according to which charged particles are characterised by a spatial distri- 
bution (spherical or other) motivated important works in classical electro-magnetism 
at the beginning of our century. Entanglement really marks a departure from such 
lines of thought, which confirms the deep intuition of Schroedinger, already men- 
tioned in the introduction, who described entanglement as the characteristic trait of 
quantum mechanics, "the one that enforces its entire departure from classical lines 
of thought" [1]. 

The present work was motivated by the results presented in the references [8,9]. 
In these papers it is argued and shown that retrievable, usable quantum information 
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can be transferred in a scheme which, in striking contrast to the quantum telepor- 
tation schemes, requires no external channel and does not involve the transfer of 
a quantum state from one subsystem to the other. Entanglement-free interaction 
between two mutually scattering particles (in the three dimensional, physical space) 
plays a crucial role in this scheme. The previous remarks suggest that localisation of 
at least one of the particles is a necessary ingredient of such protocols for quantum 
information transfer. For instance, in the test-particle limit the massive particle is 
localised while in the classical hmit, both particles are localised. It is easy to show 
that if at least one of the two interacting particles is not well localised (bilocated 
for instance), and that the particles interact through a position-dependent potential 
(action at a distance), they are highly likely to end up in an entangled state. 
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